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Abstract 


One of the dangers facing Iran’s plains is subsidence. Land subsidence due to various reasons, 
such as over-harvesting of underground water sources and climate changes, has caused many 
problems in agricultural lands, roads, and power and energy transmission lines. Therefore, it is 
important to deal with the causes and influencing factors to control the risks. Ardabil plain is also 
one of the areas prone to subsidence due to climate changes and the decrease of underground 
water in recent years. The aim of the present research is the zoning subsidence risk in this plain. 
First, the factors influencing subsidence in the Ardabil plain (slope, land use, lithology, distance 
from the fault, distance from the waterway, drop in water level, distance from the city and the 
village) were identified. Then the layers of data were applied in the geographic information 
system. In the next stage, the weighting of the investigated factors was done using ANP method 
and in the Super Decision software, and the final analysis and modeling were done using the 
MABAC method as one of the multi-criteria decision making methods. Finally, the resulting map 
was Classified into five categories from very low risk to very high risk. The results showed that 
such factors as water level drop, distance from the river and lithology have the highest weighting 
factor. Moreover, 244.29 and 370.59 square kilometers, respectively, of the area of this plain, are 
in very dangerous and dangerous classes. Finally, due to the high potential of Ardabil plain in 
terms of subsidence, protective and management measures should be taken into consideration by 
relevant authorities and institutions. 
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جغرافیا و مخاطرات محیطی. سال یازدهم. شمارة چهل و چهارم زمستان ۶۰۱ صص PW‏ 
مقاله پژوهشی 


ویژه‌نامه (چالش جهانی فرونشست زمین: مدیریت بحران يا بحران مدیریت) 


ارزيابى و بهنه‌بندی خطر فرونشست با استفاده از الگوریتم تطبیقی MABAC‏ و ANP‏ 
(مطالعة موردى: دشت اردبيل) 


(موسی عابدينى '- استاد ژئومورفولوژی» دانشكده علوم اجتماعی» دانشكاه محقق اردبيلى؛ اردبیل» ايران. 
ليلا آقايارى- دانشجوى دكترى ژئومورفولوژی. دانشكده علوم اجتماعی, دانشگاه محقق اردبیلی» اردبیل» ايران. 
صیاد اصغری سراسکانرود- استاد ژئومورفولوژی دانشکده علوم اجتماعی, دانشگاه محقق اردبیلی» اردبیل ايران. 


تاريخ دریافت: ۰۰/۱۰/۲۶ع ۱‏ تاريخ بازنگری: ۱۶۰۱/۳/۸۱ تاريخ تصویب: ۱۶۰۱/۶/۶ 


حکیده 


یکی از مخاطرات پیشروی دشت‌های کشور مخاطره فرونشست است. پدیده فرونشست زمين به دلایل 
مختلف ازجمله برداشت بیش ازحد منابع آب زیرزمینی و تغییرات جوی سبب بروز مشکلات و معضلات 
فراوان در زمین‌های كشاورزىء جاده‌هاء خطوط انتقال نیرو و انرژی شده است؛ ازاین‌رو پرداختن به علل و 
عوامل SE‏ گذار جهت کنترل و مدیریت خطر دارای اهمیت است. دشت اردبیل نيز به خاطر تغییرات 
اقلیمی و افت سطح آب‌های زیرزمینی در سال‌های اخیر یکی از مناطق مستعد جهت مخاطره فرونشست 
است. هدف تحقیق حاضرء پهنه‌بندی خطر فرونشست در اين دشت است. در اين مطالعه dal‏ عوامل موثر 
جهت ایجاد فرونشست در دشت اردبیل (شیب. کاربری اراضی, لیتولوژی. فاصله از گسل. فاصله از 


آبراهه. افت سطح آب. فاصله‌ازشهر و روستا)» شناسایی شدند و سپس نسبت به Og‏ لایه‌های اطلاعاتی در 


سامانه اطلاعات جغرافيايى اقدام كرديد. در مرحله بعد وزن‌دهی عوامل موردبررسى. با استفاده از روش 
ANP‏ و در محبط نرم‌افزار Super Decision‏ انجام گردید و تحليل و مدلسازى نهايى با استفاده از 


Email: abedini@uma.ac.ir ۰۹۱۳۵۵۵۲۲7 نویسنده مسئول:‎ ١ 
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روش MABAC‏ به‌عنوان یکی از روش‌های تصمیم‌گیری چندمعیاره. انجام شد. درنهایت. iii‏ حاصله 
در ينج رده با خطر بسیار کم تا حطر بسیار زياد طبقه‌بندی گردید. با توجه به نتایج مطالعه عوامل افت 
سطح آب. فاصله از رودخانه و لیتولوژی بیشترین ضریب وزنی را به خود احتصاص دادند. همچنین, quU‏ 
مطالعه OUS‏ داد؛ به ترتیب ۲۶/۲۹ و ۹ کیلومترمربع از مساحت اين دشت. در طبقات بسیار پرخطر 
و پرخطر قرار دارند. درنهایت می‌توان اظهار داشت. نظر به توان بالای دشت اردبیل» از لحاظ رخداد 
فرونشست. بایستی اقدامات حفاظتی» مدیریتی در سطح دشت اردبیل مورد توجه مسئولان و دستگاه‌های 
ذیربط قرار گیرد. 

کلیدواژه‌ها: مخاطرات. پدیده فرونشست. MCDM‏ دشت اردبیل. 


۱- مقدمه 

چشم‌انداز سطح زمین به‌عنوان سیستمی پیچیده حاصل تعامل فاکتورهایی نظیر فرآیندهای جغرافیایی» اقلیمی؛ 
زمانی و فعالیت‌های انسانی است. رشد کالبدی, بدون رعايت اصول برنامه‌ریزی شهری و ساخت‌وسازهایی بدون در 
نظر گرفتن پتانسیل مخاطرات موجود در پهنه‌های در معرض مخاطرات quee‏ آسیب‌پذیری و خسارت‌ها در 
محدوده‌های مورد نظر را تشدید می‌کند. عواملی که با تغيير در چگونگی ساختار و کارکرد توي وكرافى طبیعی و ريز 
اقلیم به‌ویژه در محيطهاى انسان‌ساخت. زمینه‌ساز ایجاد مخاطراتی نظیر زمین‌لغزش» سیلاب‌شهری» زمین‌لرزه و 
فرونشست می گردد (فنی و همکاران. ۱۳۹۱). براساس تعریف سازمان زمین‌شناسی ایالت متحده پدیده فرونشست 
زمين شامل فروریزش يا نشست روبه پائین سطح زمین است که می‌تواند دارای حركت قائم رو به پائین سطح زمین 
و کمی افقی باشد (عابدینی»۱۳۹۱). پدیده فرونشست زمين به دلایل مختلف ازجمله برداشت بیش از حد منابع آب 


زیرزمینی و تغییرات جوی سبب بروز مشکلات و معضلات فراوان در زمین‌های کشاورزی, جاده‌ها؛ خطوط انتقال 


سرعت شکل‌گیری» به دو دسته سریع و تدریجی تقسیم می‌گردد. حالت تدریجی به دلیل تلفات انسانی کم آن در 
مقایسه با دیگر پدیده‌های طبیعی مانند ;3« لغزش و سيلاب كمتر موردتوجه قرار فى كيرد (مدنی و همکاران. 
٦‏ ( ۱۵ 


اما فرونشست در حالت سریع باعث آسیب به سازه‌ها و تلفات انسانی می‌گردد. فرونشست زمين در سال‌های 


اخير به مشکلی جهانی تبدیل شده است. چراکه در اكثر کشورها در حال وقوع است. فرونشست سبب تخریب 
کانال‌های آبیاری» شکست خطوط لوله‌های col‏ و نفت. خطوط انتقال نیرو و تخریب سازه‌ها و تشدید فرایند زمين 
لغزه و سیلاب‌ها می‌شود (محمودی و همکاران. KITAY‏ فرايندهاى طبيعى مانند حركات تكتونيكى» فوران oltaisi‏ 


سال یازدهم ارزیابی و پهنه‌بندی خطر فرونشست با استفاده از .... £o‏ 


وسیع می‌گردند (ژویی ' و همکاران C Y‏ به‌طورکلی فرونشست متأثر از فعالیت‌های انسانی مانند برداشت بىرويه 
آب از سفره‌های آب زیرزمینی و پارامترهای زمین‌شناسی است که در صورت عدم شناخت اين پدیده سبب آسیب 
جدی به زیرساخت‌های اقتصادی می‌گردد (گالووی" و همکاران. ۲۰۰۹). پژوهشگران رابطه بين توسعه شهری و 


tij‏ زمين را رابطه‌ای دو جانبه می‌دانند که هردو بر هم تأثیر گذاشته و تأثير می‌پذیرند. به گفته آن‌ها تمدن 


شهرها به جهت رشد سریع در توسعه شهری در بخش‌های مختلف صنعت. تجارت و حمل‌ونقل اثرات 
زیست‌محیطی منفی از خود بر جا می‌گذارند؛ مانند: تغییر در کاربری اراضی کشاورزی به مسکونی» مصرف بی‌رویه 
منابع آب‌های زیرزمینی جهت گسترش فعالیت‌های صنعتی و افزايش جمعیت اين اثرات منفی درنهایت منجر به بروز 
پدیده فرونشست می‌شود و همین پدیده می‌تواند برنامه و فرایند توسعه شهری را متأثر سازد (منتظریون و همکاران؛ 
۸ 

به‌منظور شناسایی و کاهش پیامدهای ناشی از پدیده فرونشست نیاز به یک سیستم يايش فرونشست بیش‌ازپیش 
در کشور و دشت اردبیل احساس می‌شود. اولین گام در پایش فرونشست اندازه‌گیری مقدار جابه‌جایی ناشی از آن در 
سطح زمین است. با علم به میزان نرخ فرونشست و همچنین پهنه و گستره مناطق تحت تأثير, نه‌تنها قادر خواهیم بود 
میزان پیشرفت اين پدیده را در یک منطقه شناسایی کنیم و از پیشرفت بیشتر آن جلوگیری نماييم. بلکه می‌توانیم 
اطلاعات باارزشی از ویژگی‌های زمین‌شناسی و سفره آب زیرزمینی به دست آوریم. 

امروزه؛ از یک سوء سیستم اطلاعات جغرافیایی. ابزاری مؤثر برای سیستم‌های پشتیبانی تصمیم‌گیری تعاملی برای 
عملیات مدیریت ريسك هستند (بوخیر " و همکاران» (V‏ و در سال‌های اخير. جهت مطالعه فرونشست مورد 
استفاده فراوان قرار می‌گیرند. از سوی دیگر برای كشف طیف وسیعی از گزینه‌ها از نظر درگیری‌های عینی و 
معیارهای چندگانه. از روش MCDM‏ استفاده می‌شود (آهر * و همکاران» ۲۰۱۳). در اين روش‌های تصميم كيرى 
برای بهینهسازی مدلء از معیارهای مختلفی برای افزايش صحت تصمیمات استفاده می‌شود (جورجیو و همکاران" 
۵ عرب عامری و همکاران: ۲۰۱۸). لذا می‌توان گفت» استفاده از سیستم‌های اطلاعات جغرافیایی و روش‌های 
تصمیم گیری چندمعیاره با یک رویکرد تلفیقی, می‌تواند باعث تسریع در روند برنامه‌ریزی در تشخیص موارد بحرانی 
و اضطراری شده و منجر به صدور نتایج مناسبی گردد. مطالعات متعددی در نقاط مختلف جهان انجام شده و توانایی 
تحلیل‌های تصمیم‌گیری چندمعیاره» برای ارزیابی و پهنه‌بندی خطر فرونشست نشان داده است. تاکنون مطالعات 


فراوانی پیرامون پهنه‌بندی مخاطرات ازجمله فرونشست در ايران و ساير کشورها انجام شده است که به برخى از laol‏ 


1 zhu et all 

2 Galloway at all 
3 Bukheir et all 
4 Aher et all 

5 Georgiou et all 
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اشاره می‌شود. استیفن ۲ QVO‏ به بررسی علل فرونشست و عوامل تشدیدکننده و خسارات ناشی از فرونشست را 
بررسی کرده و مهم‌ترین علت فرونشست را تخلیه پیش‌ازحد آب‌های زیرزمینی و انحلال سنگ‌های کربناته اعلام کرد. 
كوبى انجوین (V‏ در کتاب fle‏ فرونشست زمین در شهر هوشی‌مین به مطالعه اين پدیده پرداخته و توصیه‌های 
مناسبی در جهت مدیریت آن در سطح شهر بیان کرده است. گامبولاتی و همکاران (۲۰۰7) به مطالعه تأثیر 
فعالیت‌های انسانی بر میزان فرونشسست زمین در مناطق مسکونی پرداختند و طبق بررسی‌های انجام شده برداشت 


بی‌رویه منابع زیرزمینی و تجمع سازه‌ها را مهم‌ترین علت فرونة = در مناطق مسکونی بیان کردند. هوانگ و 


همکاران " (۲۰۱) به بررسی فرونشست شهر دزو با روش SBAS‏ پرداختند. بررسی‌ها طی بازه زمانی دو ساله نشان 


می‌دهد که مرکز فرو بخش‌های شرقی و شمال‌شرقی با میزان فرونشست £0 میلی‌متری است و اين نشست 
زمين و افت سطح آب‌های زیرزمینی در شهرستان کرج با استفاده از روش تلفیق وزنی در محیط GIS‏ کرده و نتيجه 
کارشان نشان داده که بين وضعیت توپوگرافی. ضخامت سازند و برداشت آب ارتباط تنگاتنگی وجود دارد و بیشترین 


ميزان فرو در مناطق با برداشت obj‏ آبرفت‌های ضخیم» و مناطق يست دشت مشاهده می‌شود. عابدینی و 
همکاران (۱۳۹۵) در مورد فرونشست دشت روانسرا- سنجابی با استفاده از مدل تحلیل سلسله مراتبی کار کرده‌اند و 
نتایج نشان داده که ارتباط و همبستگی قوی ميان عمق رسوبات ریز دانه چاه‌ها و احتمال فرونشست زياد آتی وجود 
ندارد به‌طوری که حداکثر عمق رسوبات ريز دانه در نواحی مستعد فرونشینی با احتمال خطر UL‏ حدود ٩۱۲‏ متر است 
اين در حالى است که در برخی از نقاط آبخوان که استعداد کمتری برای فرونشست دارند عمق رسوبات ريز دانه بسیار 
بيشتر و حدود ۲۰ ۲۲ متر است. 

شادفر و همکاران (TAL)‏ در مقاله‌ای به پهنه‌بندی خطر فرونشست زمين با استفاده از روش AHP‏ يرداختند. 
نتایج حاصل نشان می‌دهد که فاکتور برداشت بی‌رویه با ضریب 4۸ درصد در درجه اول و عامل سنگ‌شناسی با 
ضریب VE‏ درصد در درجه دوم اهميت در ایجاد فرونشست منطقه است. همچنین نتایج نشان داد که ۵ درصد از 
مساحت منطقه در پهنه‌بندی با خطر نسبتاً زياد و خيلى زياد قرار دارند. شايان (۱۳۹۵) به پهنه‌بندی فرونشست زمین 
در حوضه آبخیز قره‌چای همدان پرداختند. جهت پهنه‌بندی از ۵ ايه ارتفاع» زمین‌شناسیء هیدرولوژیکی شبکه‌های 
زهکشی و شیب استفاده کردند. نتایج OLY‏ می‌دهد که تغییرات ناهمسان در داده‌های مورفومتری شبکه زهکشی 
مهم‌ترین فاکتور موثر در تشدید افت آب‌های زیرزمینی و وقوع فرونشست در منطقه بوده و عوامل لیتولوژیکی و 
انسانی به‌تنهایی نقشی ندارند. فنی و همکاران YAT)‏ به پهنه‌بندی و تحلیل فضایی حساسیت مخاطرات چهار گانه 


1 Stephen 

2 coein enjoein 
3 Gambolati 

4 hoang at all 
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رودخانه. کاربری اراضی و زمین‌شناسی استفاده کرده‌اند و نقشه حساسيت نهايى در قالب ۵ كلاس حساسيت 


خیلی‌زیاد» زياد متوسطء کم و خیلی کم طبقه‌بندی و به‌صورت نقشه و جدول استخراج و ارائه شده‌اند. عالی يور و 
همکاران OY)‏ به پهنه‌بندی ريسك فرونشست زمین در اثر افت سطح أب زیرزمینی با استفاده از مدل تحلیل 
سلسله مراتبی فازی دشت اردبیل پرداختند و اقدام به تهیه نقشه احتمال فرونشست زمین و تعیین مناطق دارای 
ريسك پرداختند. لایه‌های اطلاعاتی شامل ٩‏ لایه زمین‌شناسی. فاصله از گسل, افت سطح‌آب زیرزمینی» شيب و 
ضخامت آبرفت به‌عنوان معیارهای مؤثر مورد استفاده قرار گرفتند و نتایج نشان داد كه ۱۳ درصد منطقه دارای ريسك 
خيلىبالا است و حدود 1۷/۵ درصد مناطق جنوبی در بخش ریسک‌های خیلی بالای فرونشست. به دلیل استفاده 
بیش‌ازحد از آب‌های زیرزمینی است. بهیاری و همکاران (YS‏ به ارزیابی تأثیر ساختارهای فعال بر خطر 
فرونشست زمين با استفاده از مدل‌های تصمیم‌گیری چندمتغیره پرداختند. در اين مطالعه از روش AHP‏ به‌منظور 
تجزیه‌وتحلیل خطر فرونشست در منطقه موردمطالعه با استفاده از لایه‌های مطالعاتی ارتفاع, شیب. آبراهه‌ها و 


زمین‌شناسی, نقشه پهنه‌بندی das‏ فرونشست در منطقه تهیه گردید. اين نقشه OUS‏ می‌دهد كه خطر فرونشست در 
بخش جنوبی منطقه نسبت به ساير مناطق بیشتر است. كوه بنانی و همکاران (۱۳۹۸) به پهنه‌بندی گستره خطر 
فرونشست زمین با بهره‌گیری از تداحل‌سنجی راداری در دشت کاشمر و خلیل‌آباد با استفاده از تصاویر سنتینل ۱ 
a byy‏ مال هاف aul dosis AY UAE‏ نشان داف این dice‏ جر pagi ca ۱۷ Gib cols‏ در سنال 


فرونشست دارد و پیشتر در حومه جنوب غرب شهر کاشمر و مناطق غربی دشت خليل أباد متمرکز است. منتظریون و 


دشت آبخوان نورآباد با استفاده از روش تداخل‌سنجی 
راداری پرداختند و جهت بررسی وضعیت افت سطح آب‌های زیرزمینی از روش درونیابی IDW‏ استفاده کرده‌اند و 
نقشه‌های حاصل از فرونشست نشان می‌دهد که نرخ نشست در هرسال به‌طور میانگین ‏ سانتی‌متر و طی بازه زمانی 
٤‏ ساله ۱١‏ سانتی‌متر است. ابراهیمی و همکاران (۱۳۹۹» به يايش میزان فرونشست محدوده شهری دشت پاکدشت 
با روش تداخل‌سنجی پرداختند. نتایج حاصل از روش تداحل‌سنجی راداری محدوده مطالعاتی در طی بازه زمانی 
یک‌ساله بين ۱۵ تا ۸۵ میلی‌متر فرونشست را OUS‏ می‌دهد و مهم‌ترین و موثرترین عامل در فرونشست را هم افت 
سطح آب‌های زیرزمینی عنوان کرده‌اند. شاه کرمی و همکاران (۱۳۹۹) مناطق فرونشست در اثر اضافه برداشت آب 
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زیرزمینی در دشت اراک را موردبررسی قرار دادند. نتایج به دست آمده Sly‏ اين موضوع هستند که نواحی مرزی و 
حاشیه‌ای دشت اراک در معرض افت بیش‌تر آب قرار دارند و تمرکز رسوبات ریزدانه رسی در قسمت‌های غربی 
به‌مراتب بیشتر از نواحی دیگر است. درنتیجه احتمال می‌رود اين نواحی در معرض خطر فرونشست بیشتری باشند. 


زارعی و همکاران (۱۳۹۹) به تعيين رابطه فرو 


زمین و افت سطح آب زیرزمینی با دو روش تداخل‌سنجی 
راداری و ایستگاه سایت gps‏ دشت سلماس پرداختند و نتایج OUS‏ داد که به ازای یک متر افت سطح آب زیرزمینی؛ 
۸ متر سطح زمین در دشت نشست می‌کند. جوادی Or)‏ در رساله خود به برآورد و تحلیل میزان فرونشست 
دشت همدان- بهار با استفاده از تصاوير راداری و سری زمانی SBAS‏ در طی بازه زمانی ۵ ساله و میزان فرونشست 
اين دشت را ۵۱ میلی‌متر به‌صورت سالانه برآورد کرده است و بیش‌ترین میزان فرونشست مربوط به مناطق دارای 


افت منابع آب زیرزمینی است. فرزین کیا و همکاران ( 6۱۶۰۰ به پهنه‌بندای خطر فرونشست زمین در دشت جوین با 


استفاده از مدل تحلیل شبکه‌ای فازی پرداختند در اين تحقیق از روش ANP‏ و فازی. à‏ 25 ت به‌صورت 
کلاس‌بندی برای نشان دادن ویژگی‌های مناطق بحرانی فرونشست و اولویت‌بندی مناطق مختلف از نظر استعداد برای 
وقوع پدیده فرونشست به دست آمد و نتایج بیانگر اين است که OV‏ درصد قسمت مرکزی دشت جوین با خطر زياد 
تا خیلی زياد بهحساب می‌آیند. 


هرچند که تحقیقات معدودی در مورد دشت اردبیلی انجام شده ولی موضوع آن‌ها متفاوت بوده اسست. 


YA) e‏ در مورد ple‏ فرونة 2 دشت اردبیل کار 02,5 و مهم‌ترین fle‏ أن را برداشت بی‌رویه آب‌های 
زیرزمینی» فعالیت‌های تکتونیک حاصل از گسل‌های حواشى دشت و سنگینی ناشى از انباشت آبرفت‌ها عنوان كرده 
است 
سسا 


با توجه به منابع بررسی شده محققان مختلف اين تکنیک را بهعنوان ابزار نيرومندى برای مطالعه پدیده‌های 
مختلف زمین‌شناسی همچون فرونشست مطرح كرده است. روش MABAC‏ كه براساس رويكرد برنامه‌ریزی 
مطلوب و بهینه‌سازی چندمعیاره پایه‌ریزی شده است. به‌عنوان روش تصمیم‌گیری چند معیاره در تلفیق با فرایند 
تحلیل شبکه به‌منظور پهنه‌بندی پتانسیل وقوع فرونشست دشت اردبیل استفاده شد. طبق بررسی als‏ منتشر شده 


تاکنون در دشت اردبیل با اين روش و داده‌های ورودی» پژوهشی ارائه نشده است لذا جهت برآورد دقیق‌تر میزان 


فرو و عوامل تأثیرگذار بر آن لزوم و خلأ انجام تحقیق با این روش كاملاً احساس می‌شود. لذا اهداف تحقیق 
حاضر بررسی فرونشست دشت اردبیل با استفاده از تکنیک تداخل سنجی راداری, پهنه‌بندی فرونشست. اندازه‌گیری 


معیارهای مؤثر بر فرونشست و بررسی مژثرترین عوامل در فرونشست دشت اردبیل است. 
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۲- مواد و روش‌ها 
۱-۲- منطقه موردمطالعه 


اردبیل واقع شده است. در بررسی زمین‌شناسی دشت اردبیل با مختصات جغرافیایی EA‏ درجه تا ÉA‏ درجه و ۳۰ 


دقيقه طول خاوری و ۳۸ درجه و ۲۰ دقيقه شمالی قرار دارد. زمین‌های حاشیه دشت از اطراف به سمت دشت شیب 
دارند. ارتفاع متوسط اراضی کشاورزی حدود ۱۳۲۰ متر از سطح دریاهای آزاد قرار دارد. وسعت آن حدود ۸۲۰ 
کیلومترمربع است و جزئی از حوضه آبخیز رودخانه قره‌سو به شمار می‌رود. دشت اردبیل متشکل از ۸۸ واحد 
روستایی ده واحد دهستان و ۲ واحد شهری اردبیل و آبی بیگلو است. شهر اردبیل مرکز استان اردبیل در غرب اين 
دشت ally‏ شده است. در محدوده دشت اردبیل تعداد ۲۲۵۳ حلقه چاه عمیق و نیمه‌عمیق, ۲۰ دهنه چشمه و VA‏ 


رشته قنات وجود دارد كه به مصارف کشاورزی اختصاص دارند (رستم‌زاده و ATE iS‏ 
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شکل ۱- نقشه موقعیت دشت اردبیل 
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-Y-Y‏ داده‌ها و ابزارهاى مورد استفاده 

در راستای انجام اين پژوهش از داده و تصاویر ماهواره‌ای و نرم‌افزارهای متعددی استفاده شده است که مشتمل بر 
موارد زیر خواهد بود: 

پژوهش حاضر از نوع کاربردی بوده و روش تحقیق آن تحلیلی مبتنی بر آنالیز داده‌ه؛ سیستم اطلاعات جغرافیایی و 
استفاده از فنون تحلیل چندمعیاره اسست. جهت پردازش تصاویر و تجزیه‌وتحلیل داده‌ها نیز از نرم‌افزارهای 
Excel ; Arc GIS Ecognition ENVI‏ بهره‌گیری شده است. جهت بررسی خطر فرونشست. در مرحله جمع‌آوری 
دادههاء ابتدا بايد ملفه‌های تأثیرگذار مشخص شده و بر مبنای laol‏ داده‌ها و اطلاعات مورد نیاز گردآوری و طبقه‌بندی 
شوند. بدین منظورء در اين مطالعه؛ ابتدا عوامل مؤثر (شامل: شیب لیتولوژی, کاربری اراضی فاصله از آبراهه» فاصله 
از كسلء افت سطحآب و فاصله از شهر و روستا) با توجه به شرایط طبیعی و انسانی منطقه شناسایی گردید. در 
مرحله بعد لایه‌های اطلاعاتی مربوط به هر یک از عوامل, در محیط سيستم اطلاعات جغرافیایی تهیه گردید. LU‏ 
اطلاعاتی شبکه آبراهه با استفاده از نقشة رودخانه‌های استان اردبیل» تهیه گردید. لایه‌های شیب نیز با استفاده از مدل 
رقومی ارتفاعی, دانلود شده از سايت (USGS)‏ تهیه شد. لایه‌های اطلاعاتی مربوط به لیتولوژی (مقاومت سنگ‌ها) و 
كسلهاء با رقومی‌سازی از روی نقشه زمین‌شناسی اردبیل؛ با مقیاس ۱:۱۰۰۰۰۰ تهیه شد. برای استخراج کاربری 
اراضی محدوده موردمطالعه» از تصاویر سنجنده لندست هشت (سنجنده (OLI‏ به شماره گذر VW‏ و ردیف ۳۶ 
تاريخ ۰ ۷ dl Y‏ شده از «(earthexplorer.usgs.gov ) cole‏ استفاده شده است (حدول ۱). 

تصاوير ماهوارهاى ايتيكال: تصوير ماهوارهاى ايتيكال در اين تحقيق بهمنظور استخراج نقشه كاربرى اراضى سال 
۰ استفاده شد. مشخصات تصاوير لندست اخذ شده از سايت <(earthexplorer.usgs.gov)‏ در جدول ١‏ ارائه 
شده است. به‌منظور تصحيحات اتمسفرى و راديومتريكى از نرم‌افزار ENVIS.3‏ و به‌منظور استخراج نقشه‌های مربوطه 
از نرم‌افزار ARCGIS10.5‏ استفاده شد. برای طبقه‌بندی کاربری اراضی از روش طبقه‌بندی شیءگرا در 


نرم‌افزار developer64‏ 0 استفاده شد. همچنین به‌منظور به دست آوردن صحت طبقه‌بندی از Google‏ 


Earth‏ استفاده شد. 
حدول ۱- اطلاعات تصاویر سنحنده لندست ۸ 
SCENE-CENTER-TIME SUN-AZIMUTH ROW PATH Sansor name Year‏ 
VT \£0,0+ ۷ yt vw Landsat 8 ۳۹۰‏ 


در روش طبقه‌بندی شیءگرا؛ اطلاعات طیفی با اطلاعات مکانی ادعام گردیده و پیکسل‌ها پراساس شکل. بافت و 
تن خاکستری در سطح تصوير با مقیاس مشخص, قطعه‌بندی شده و طبقه‌بندی تصوير براساس اين قطعه‌ها انجام 
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می‌شود. در قطعه‌بندی پیکسل‌ها توسط الگوریتم‌های مختلف. با نسبت‌های مختلف طیفی و شکلی طبقه‌بندی و 
براساس ویژگی‌های ab‏ و فضایی در قالب اشيا گوناگون دسته‌بندی شده‌اند. طی اين فرآیند شی‌های تصویری 
متناسب با همگن بودن يا ناهمگن بودن براساس پارامترهای مقیاس» رنگ» شکل» ضریب نرمی و شکل فشردگی 
اناد هد gray)‏ لا ۱۳۹۶ با آرمون y‏ عطا بهترین انداده که تان steal Lob dta S‏ مخقلت go‏ باد 
تعیین شد. در اين پژوهش تصوير قطعه‌بندی شده دارای مقیاس LO‏ شکل ۰/۵ ضریب فشردگی ۰/۵ انجام گردید. 
بعد از مشخص كردن تعداد و کلاس‌های موردنظر در طبقه‌بندی تعيين شی‌های نمونه‌های آموزشی انجام و در 
نرم‌افزار بر سطح تصوير اعمال گردید. سپس طبقه‌بندی با استفاده از نزدیک‌ترین همسایه و روش درخت 
تصميم كيرى صورت گرفت. لازمه استفاده از هر نوع اطلاعات موضوعی, آگاهی از صحت و درستی OT‏ است. دقت 
ouo GI dls sees‏ کی A‏ هل و شیب jd EA‏ 

فرآیند تحلیل شبکه‌ای (ANP)‏ 

روش (ANP)‏ شکل توسعه یافته‌ای از روش (AHP)‏ است که قادر است همبستگی و بازخوردهای موجود بين 
pole‏ در يك تصمیم‌گیری را مدل‌سازی نموده و تمامی تأثیرات درونی اجزای موثر در تصمیم‌گیری را منظور و وارد 
محاسبات نماید. در مدل AHP‏ تنها معیارها با معیار بالایی خود ارتباط داشتند؛ یعنی معیارها به هدفهاء زیرمعیارها 
به معیارها و گزینه‌ها هم به زیرمعیارها وابسته بودند» حتی خود معیارها هم به همدیگر وابسته نبودند. ولی در مدل 
ANP‏ نهتنها خوشه‌ها بر عناصر و عناصر بر گزینه. گزینه بر عناصر تأثیر می‌گذارند. Sh‏ حتی pole‏ بر خودشان و 
بر So‏ حوشه‌ها نیز اثرگذارند. به‌طو ركلى» روش (ANP)‏ دارای مراحلی نیز هست که شامل. ۱- تعيين معیارها و 
شاخص‌ها ۲- تعیین روابط و ارتباطات بين عناصر و خوشه‌ها: در اين گام» es‏ تصميمكيرى به ساختار شبکه‌ای 
تجزیه می‌گردد.۳- مقایسات زوجی بين عناصر و خوشه‌ها: ماتریس‌های مقایسات زوجی مشابه با روش (AHP)‏ 
است و زوج‌های عناصر تصمیم‌گیری در هر خوشه به نسبت اهمیتشان در جهت شرطهای JES‏ آن‌ها مقایسه 
می‌شوند. 6- تشکیل سوپرماتریس‌ها: الف- تشکیل سوپرماتریس اولیه يا غير وزنی ب- تشکیل سوپرماتریس وزنی 
ج- تشکیل سوپرماتریس حدی: سوپرماتریس وزنی» به توان حدی می‌رسد تا عناصر ماتریس همگرا شده و مقادیر 
سطری bol‏ هم برابر شوند. ۵- نتایج ماتریس خوشهها. T‏ نتیجۂ نهایی عناصر: ضرایب سوپرماتریس در ضرایب 
ماتریس خوشه‌ها نرمال شده و درنهایت. نتيجة نهایی عناصر و اولویت آن‌ها مشخص می‌گردد (فرجی سبكبار و 
همکاران ۱۳۸۹). 

گام‌های روش ماباک (MABAC)‏ 


۱- تعيين معیارها و گزینه‌های پژوهش. Y‏ تشکیل ماتریس تصمیم: ماتریس تصمیم در اين روش به‌صورت 
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سلول نيز در واقع امتیاز هر گزینه نسبت به هر معیار است. ۳- در اين مرحله نسبت به استاندارسازی و يا 
به‌عبارت‌دیگر بی‌مقیاس نمودن معیارهاء اقدام می‌گردد. در اين پژوهش, از روش فازی و تابع عضویت خطی 
(Linear)‏ و سیگموئیدال (Sigmoidal)‏ به‌منظور استانداردسازی استفاده شده است. در مجموعه‌های فازی, مقدار 
يك به‌عنوان بیشترین ارزش به حداکثر عضویت و صفر به‌عنوان کمترین ارزش» به حداقل عضویت لحاظ می‌گردد 
(حسینی و همکاران (MA‏ 

6- وزن‌دار كردن ماتریس نرمال: در اين گام با استفاده از رابطه زیر ماتریس نرمال را وزن‌دار مىكنيم (رابطه AY‏ 
در مطالعه حاضر جهت وزن‌دهی معیارها از روش ANP‏ استفاده شده است و نتایج حاصل از وزن‌دهی ANP‏ با 
توجه به رابطه زین به تک تک لایه‌های اطلاعاتی لحاظ گردیده است. 

(V_ij=W_j*(nij + 1 ) ۱ رابطه‎ 

که در آن nij‏ عناصر ماتریس نرمال شده N)‏ را Wiss p OUS‏ ضرایب وزنی معیار را OUS‏ می‌دهد. 

۵- تعيين مرز ناحيه شباهت ماتریس (ع): در رابطه vijé‏ عناصر ماتریس وزنی (V)‏ و m‏ تعداد کل گزینه‌های 
جایگزین را نشان می‌دهد. پس از محاسبه مقدار gi‏ ماتریس محدوده‌های تقریبی مرزی با توجه به معیارهای LG‏ 
فرمت ۸*1 شکل م ىكيرد و 7 تعداد کل معیارهایی را که گزینه‌های ارائه شده برای آن‌ها انتخاب می‌شود را OUS‏ 
می‌دهد. 

رابطه G-D^mv ij (Vm) Y‏ ]120 ع) 

بنابراین اگر 0 معیار داشته باشید یک ماتریس Glxn‏ به‌صورت زیر eal‏ داشت (رابطه ۳): 


[G-Igl.g2.....gn Y رابطه‎ 

aal uale‏ گر ها قاس اھک ور ای شش با o las ala: ada Sana‏ اسرد كيه 
دست می‌آید. در واقع بايد ماتریس وزن‌دار را از ماتریس 8 کم کرد. 

Q-V-G ٤ رابطه‎ 

(GI وه بای ات‎ Laos, اوو بای‎ Ale sl 1: قواق‎ Oo e Dae sess Jos 
وضعیت هر گزینه را مشخص کرد.‎ 

۷- انتخاب گزینه بهینه: در روش ماباک MABAC‏ با استفاده از رابطه ف امتیاز نهایی هر گزینه را مشخص کرده 
تس کرت وای ی ر 

S_FZ(qij);=1.2....n0j=1.2....m ۵ رابطه‎ 


سال یازدهم ارزیابی و پهنه‌بندی خطر فرونشست با استفاده از .... oy‏ 


محاسبه مقادیر توابع kaa‏ توسط گزینه‌ها به‌عنوان مجموع فاصله‌های جایگزین از مناطق تقریبی مرزی 41 به 
رابطه n‏ تعداد معیارها را نشان می‌دهد و m‏ تعداد گزینه است (علیزاده و همکاران (Ye)‏ 


جدول ۲- نوع توابع فازی‌سازی لایه‌های موضوعی تأثیرگذار بر خطر فرونشست دشت اردبیل 


متغیر توضیحات نوع alt‏ فازی 
كيت کلاسه‌بندی مجدد و اختصاص کدها بر اساس اهمیت آن در وقوع خطر فرونشست و سپس فازی کردن. | تابع‌بزرگ فازی 
کاربری اراضى اختصاص كدها بر اساس اهميت آن در وقوع خطر فرونشست و سپس فازى كردن. تابع‌بزرگ فازى 
لیتولوژی اختصاص كدها بر اساس اهميت آن در وقوع خطر فرونشست و سپس فازى كردن. تابع بزرگ فازى 
فاصله از asl yl‏ در مناطق مجاور رودخانه. پتانسیل وقوع مخاطره فرونشست کاهش می‌یابد. تابع نزدیک فازی 
فاصله از گسل در مناطق مجاور گسل. پتانسیل وقوع مخاطره فرونشست افزایش می‌یابد. تابع نزدیک فازی 
فاصله از شهر و 
5 در مناطق مجاور شهر و روستاء پتانسیل وقوع مخاطره فرونشست افزايش می‌یابد. تابع نزدیک فازی 
زوس 
افت سطح آب با کاهش میزان سطح آب جاههاء خطر فرونشست افزایش مى يابد. ati‏ خطى کاهشی 


معیارهای مطرح در پهنه‌بندی فرونشست دشت اردبیل 

لیتولوژی: قدیمی‌ترین سازندهای شناخته شده در دشت اردبیل دولومیت و دولومیت Sal‏ متعلق به ژوراسیک 
است. رخساره‌های کرتاسه ماسه سنگ و کنگلومرای قرمز می‌باشد. رخساره‌های دوران سوم شامل ماسه‌سنگ‌های 
توفی و توفی برشى با ميان لایه‌های گدازه آندزیتی» آهک كريستاليزهه آندزیت بازالت با پیروکسن زياد مربوط به 
پالئوژن می‌باشد. آبرفت‌های کواترنر که شامل تراس‌های قدیمی» رسوبات جوان رودخانه‌ای» رسوبات آهکی تراورتن» 
توف» lal ASE‏ مارن و رسوبات جوان دشت بیشترین سهم را در دشت اردبیل و بخش مرکزی أن که بیش‌ترین 
ميزان تمرکز چاه‌ها و قنات‌ها را دارده از رسوبات آبرفتی تشکیل شده است كه در مرکز و قسمت‌شرقی دشت 
بيشترين میزان فرونشست دیده می‌شود و در ساير بخش‌ها به‌واسطه شیب بیش‌تر و نبود رسوبات آبرفتی زمینه برای 
وقوع فرونشست بسیار کمتر از ساير بخش‌ها است (شکل ۲ الف». 

شیب: دشت اردبیل در شمال‌غرب ايران و در قسمت‌شرقی فلات آذربایجان به‌صورت حوضه بسته‌ای است که از 
شمال و غرب منطقه به كوه سبلان و کوه‌های اطراف أن و از شرق و جنوب به رشته‌کوه‌های تالش منتهی می‌گردد؛ 
بنابراین زمین‌های حاشیه دشت از اطراف به سمت دشت شیب دارند. بررسی وضعیت ارتفاع OUS‏ می‌دهد که در 


بخش مرکزی دشت کم‌ترین ميزان ارتفاع وجود دارد. شواهد فرو نيز بیش‌تر در اين طبقه ارتفاعی دیده 


می‌شود. در مورد لايه شیب. محدوده شیب ۰-۵ کمترین میزان شيب در محدوده منطقه می‌باشد که اين شيب نقش 


oí‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


مؤثرى در جمع‌آوری منابع آب زیرزمینی فراهم می‌آورد. هر جه میزان شيب افزایش می‌یابد. از میزان تمرکز چاه‌ها 
کاسته و میزان فرونشست زمین هم GRAS‏ می‌یابد (شکل Y‏ ب). 

فاصله از رودخانه: عامل نزدیک به رودخانه به‌عنوان عامل مزاحم در ایجاد فرونشست شناسایی گردید (شایان و 
همکاران. ۱۳۹۵) و با فاصله گرفتن از رود امتیاز بیشتری به زمین‌های منطقه تعلق م ىكيرد. هر جه از آبراهه‌ها فاصله 
می‌گیریم نیاز امستفاده از آب‌های زیرزمینی به دلیل عدم وجود و استفاده از آب‌های سسطحی بالا می‌رود و در پی 
استفاده نادرست از آب‌های زیرزمینی و افت سطح آب‌های زیرزمینی» میزان احتمال فرونشست افزایش مى يابد 
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شکل ۲- نقشه معیارهای مطرح در پهنه‌بندی فرونشست دشت اردبیل؛ الف) نقشه لیتولوژی. ب) نقشه شیب. ج) 


نقشه آبراهه 


فاصله از شهر و روستا: عوامل انسان‌ساحت مانند زمین‌های کشاورزی ابنیه بزرگ و صنعتی» تأسیسات زیرزمینی 
مانند تونل» خطوط انتقال نیرو و عملیات حفاری و معدن کاری» سکونت‌گاه‌های شهری و روستایی و شبکه‌های 
حمل‌ونقل جاده‌ای با افزايش مصرف آب‌های سطحیی زیرزمینی» تحمیل بار وارده بر زمين در اثر احداث 
ساختمان‌های بزرگ و تراکم زیرساخت‌ها میزان آسیب‌پذیری اين واحد را افزايش داده و سبب تشدید فرونشست 
شده‌اند (شکل ۳ الف). 

بررسی تراکم چاه‌ها و افت سطح آب: براساس نتایج ارائه شده از مطالعات محققان, به نظر می‌رسد که اعمال 


فشار بي شازحد بر آبخوان ممکن است موجب بروز پدیدۀ فرونشست شود. ازاین‌رو به بررسی وضعیت سفره آب 


سال يازدهم ارزیابی و پهنه‌بندی خطر فرونشست با استفاده از .... 00 


زيرزمينى در محدودهٌ موردمطالعه پرداخته شد. بهمنظور بررسى و ارزيابى نوسانات سطح آب زیرزمینی» از آمار 
چاه‌های بيزومترى موجود در دشست. ميزان افت سطح آب زيرزمينى طى سالهاى WV‏ تا ۱۳۹۹ محاسبه شد. طبق 
بررسی‌های انجام شده مشخص است که سطح col‏ دشت. رو به پایین رفتن است. به‌طوری‌که کمترین افت سطح 
آب ۱/۳ متر و بيشترين ۲۳ متر می‌باشد. علت اين افت» برداشت‌های غیراصولی به دلیل افزايش حفر چاه‌های عمیق 
و نیمه‌عمیق به‌صورت مجاز و غیرمجاز: کاهش بارندگی در سال‌های اخير و به‌تبع آن کاهش آب‌های سطحی در 
محدوده موردمطالعه می‌باشد. براین‌اساس بیش‌ترین افت در بخش‌های مرکزی و شرقی دشت مشاهده می‌شود که 
احتمال وقوع فرونشست در این قسمت از دشت افزايش مىيابد (شکل ۲ ب). 

فاصله از گسل: به دلیل اثر گذاشتن فعالیت گسل‌ها (نظیر زلزله» بر روی فرونشست زمین» بررسی گسل‌های 
منطقه ضرورت می‌یابد. گسل‌ها به نوبه خود نقش مؤثر و مهمی در ميزان توسعه فرونشست در سطح دشت دارد. 
به‌گونه‌ای که هر جه میزان فاصله از حطوط گسل بیش تر باشد نشان‌دهنده اين است که منطقه تناسب کمتری برای 
احتمال فرونشست زمين را دارا بوده و هر جه میزان فاصله به خطوط گسل نزدیک‌تر شود احتمال وقوع فرونشست 
بیش‌تر می‌شود. در دشت اردبیل نيز به علت وجود گسل‌های فراوانی نظیر اردبیل» بزغوش, آستارا و.؛ که دارای روند 


شمال خاورى - جنوب باختری می‌باشد و به همین دلیل ضریب آسیب‌پذیری بالایی دارد (شکل ۳ ج). 
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شکل -Y‏ نقشه معیارهای مطرح در پهنه‌بندی فرونشست دشت اردبیل؛ الف) نقشه شهر و روستا ب) نقشه افت 


o‏ جغرافیا و مخاطرات محیطی شمارة چهارم 


کاربری اراضی: انواع کاربری در دشت اردبیل به‌صورت نواحی انسان‌ساخت» مرتع خوب. مرتع متراکم مرتع 
متوسط مرتع ضعیف پوشش برف دائمی. دیم‌زار. مناطق کشاورزی» زمين بای باغات و پوشش جنگلی» شهرک 
صنعتی» مناطق آبی و اراضی pl‏ می‌باشد. ازآن‌جهت دشت اردبیل دارای توپوگرافی هموار می‌باشد و از آپرفت‌های 
جوان و تراس‌های قدیم و جدید در اثر ته‌نشست آبرفت‌های رودخانه‌ای با شیب کم تشکیل شده است محل مناسبی 
برای فعالیت‌های زراعی و باغات ایجاد کرده است. لذا اغلب اراضی کشاورزی و نيز باغات در اين Joly‏ گسترده 
شده است. بررسی معیار کاربری اراضی SE‏ از اين موضوع است که کاربری اراضی و Sle‏ و نواحی انسان‌ساز و 
شهرک صنعتی بيش ترين درصد از مساحت مناطق پرخطر را به خود اختصاص داده است؛ که با مصرف بالای منابع 


آب زیرزمینی جهت مصارف روزانه و آبیاری احتمال وقوع فرونشست را افزایش می‌یابد. 
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شکل ۶- نقشه معیارهای مطرح در پهنه‌بندی فرونشست دشت اردبیل؛ نقشه کاربری اراضی 


۳- نتایج و بحث 


ازاین‌ری لازم است معیارها قبل از ترکیب با یکدیگر استاندارد شود. در اين مطالعه. مرحله مربوط به ارزش گذاری و 


استانداردسازی به‌صورت توأم و بر مبنای ارزش عضویت در مجموعه فازی در نظر گرفته شده است. در پژوهش 


سال يازدهم ارزیابی و پهنه‌بندی خطر فرونشست با استفاده از .... ov‏ 


حاضس ارزش عضویت بين مقادیر ١ Ge‏ قرار دارد. نقشه‌های استاندارد de‏ مربوط به معيارهاى مطرح در مطالعه 
حاضر در شکل ۵ نشان داده شده است. 
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aT‏ لاريريى اراضى نفشه اسداندارد شده اسل از شهر و روسنا 
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شکل ۵- نقشه‌های استاندارد شدۀ فازی معیارهای مطرح در پهنه‌بندی فرونشست دشت اردبیل 


0۸ جغرافیا و مخاطرات محیطی ob‏ چهارم 


غیروزنی يا اولیه مربوط به مدل» سوپرماتریس غیروزنی و سوپرماتریس حدی» آخرین مرحله برای تعيين ارزش و 
ضریب نهایی عناصرء محاسبة نتایج ماتریس خوشه‌ها و نرمال‌سازی ضریب pole‏ در سوپرماتریس حدی توسط 
ضریب خوشه‌ها می‌باشد. نتيجة مقایسات زوجی بين گروه‌ها که حدول اوزان گروه‌ها نامیده می‌شود. در حدول )« 


آورده نله Cael‏ 


جدول ۳- ماتریس اوزان گروه‌ها 


خوشه‌ها انسانی زمین‌ساختی فیزیوگرافی هیدرولوژی 

۳۲ vvv Zn ۰,۲ انسانی‎ 
SP D OM SPA زمین‌ساختی‎ 
JAM efe AN JANE فیزیوگرافی‎ 
۳۲ ۳۳۳ Ur 0 هیدرولوژی‎ 


نتيجة نهایی عناصر: ضرایب سوپرماتریس در ضرایب ماتریس خوشه‌ها نرمال شده و درنهایت. نتيجة نهایی 
عناصر و اولویت آن‌ها مشخص شد (جدول (E‏ كه معیارهای افت سطح آب چاه‌های موجود در منطقه» فاصله از 


رودخانه» لیتولوزی و کاربری اراضی به ترتیب بیشترین تأثیر را در فرونشست دشت اردبیل را داشتهاند. 


حدول 5- وزن (alg‏ حاصل از روش ANP‏ 


معیار فاصله از شهر و روستا | کاربری اراضی فاصله از گسل لیتولوژی شیب | افت سطح آب | فاصله از رودخانه 


وزن نهایی ۷ 4۷ (Y‏ 10۲ ۰/۱ ۰/۲۳۱ 2۳۹ 


پس از اعمال كردن وزن‌های حاصل از روش ANP‏ در هر يك از نقشه‌های استاندارد شده و پیاده‌سازی دیگر 
مراحل مدل MABAC‏ نقشه نهایی پهنه‌بندی فرونشست در ۵ Ab‏ بسیار پرخطر تا بسیار کم cda‏ به دست آمد 
(شکل ۱). با توجه به نتایج به دست آمده 6 كيلومترمربع از مساحت شهرستان در Abb‏ بسیار پرخطر و 
۵۹ کیلومترمربع. در طبقه پرخطر قرار دارد. همچنین ۳۸۳/۷۵ کیلومترمربع از مساحت دشت. دارای خطر متوسط 
و به ترتیب 14۱/۲۳ و ۲۸۱/۷۱ کیلومترمربع دارای خطر بسیار کم و کم می‌باشند (جدول ۵). 
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شکل - نقشه پهنه‌بندی خطر فرونشست دشت اردبیل 


جدول ۵- مساحت مربوط به هر یک از طبقات خطر 


طبقه خطر بسیار پرخطر پرخطر خطرمتوسط كم خطر بسیا ركم خطر 
مساحت به کیلومترمریع ۱۹ ۳۷/۰/۵۹ ۳۸۵ 11/۲۹ AW‏ 
مسالحت به درد vv Yo/ov YYAA 1/4۷ VEA‏ 


طبق بررسىهاى انجام شده و نتايج به دست آمده ۳۵/۵۵ درصد از منطقه موردمطالعه در يهنه بسيار يرخطر و 
يرخطر قرار دارد. با بررسى و محاسبه موردى معيارهاى تأثي ركذار بر فرونشست و يراكنش آن‌ها در يهنه بسيار يرخطر 
می‌توان بیان كرد كه فاصله از رودخانه كه طبق جدول (E)‏ مهمترين و تأثیرگذارترین معيار در فرونشست منطقه 
موردمطالعه مىباشد و هر جه از آبراهه‌ها فاصله می‌گيريم نياز استفاده از آب‌های زيرزمينى به دليل عدم‌وجود و 
استفاده از آب‌های سطحى بالا مىرود و در بى استفاده نادرست از آب‌های زيرزمينى و افت سطح آب‌هاء ميزان 
احتمال فرونشست افزایش می‌پاید. لذا ۶۲/۳۲ کیلومتر از منطقه بسیار پرخطر در فاصله ۱۸۰۰۰ متری از cota p;‏ 
۸ کبومتر در فاصله ۰ متری و ۱۱/۷۱ کیلومترمربع در فاصله ۰ مترى از رودخانه قرار دارد. هر جه 
میزان فاصله به حطوط كسلى نزدیک‌تر باشد احتمال وقوع فرونشست افزايش می‌یابد و هر جه میزان فاصله از 
خطوط گسلی بیشتر باشد نشان‌دهنده اين است که منطقه تناسب کم‌تری را برای احتمال فرونشست دارا می‌باشد. 


1 جغرافيا و مخاطرات محیطی شمارة چهارم 


۱ کیلومتر از پهنه بسیار پرخطر در فاصله ۵۰۰۰ متری از گسل» ۸۰/۹۸ کیلومتر در فاصله ۱۰۰۰۰ متری»۲۱/۳۵ 
کیلومتر در فاصله ۱۵۰۰۰ متری و ۰/۱ کیلومتر در فاصله ۲۰۰۰۰ متری از گسل قرار دارد. بيشترين احتمال وقوع 
فرونشست در محدوده شیب بين ۰ تا ۵ درصد اتفاق می‌افتد زيرا اين شیب نقش مؤثرى در جمع‌آوری منابع آب 
زیرزمینی فراهم می‌آورد. هر جه میزان شیب افزایش می‌یابد از میزان تمرکز چاه‌ها کاسته و میزان فرونشست زمین هم 
کاهش مىيابد. ۱۲۷ کیلومتر از منطقه بسیار پرخطر در محدوده شیب ٠‏ تا ۱۵ درصدء ۸۱/۹ کیلومتر در شیب ۱۵ تا 
۰ درصد و ۳۵/۱ کیلومتر در محدوده شیب ۳۰ تا £0 درصد قرار دارد. عوامل انسان‌ساخت مانند ابنیه‌های بزرگ و 
صنعتی. تأسیسات زیرزمینی» عملیات حفاری, شبکه‌های حمل ونقل جاده‌ای و سکونت‌گاه‌های شهری و روستایی با 
تحمیل بار وارده بر زمین میزان آسیب‌پذیری را افزايش داده و سبب تشدید فرونشسست می‌شود به همین دليل 
بيش ترين میزان فرونشست در پهنه بسیار پرخطر در فاصله ۲۰۰۰۰ متری از شهر و روستا با ۱۳۶/۰۹ کیلومتر در 
فاصله 2۰۰۰۰ متری از ۸۰/۲۰ کیلومتن در فاصله ۰۰۰۰ متری ۲۱ کیلومتر و در فاصله ۸۰۰۰۰ متری 51 کیلومتر 
قرار دارد. آبرفت‌های کواترنر كه شامل تراس‌های قدیمی. رسوبات جوان رودخانه‌ای» رسوبات Sal‏ تراورتن» توف 
خاكسترء لاهار, مارن و رسوبات جوان دشت بیشترین سهم را در دشت اردبیل و بخش مرکزی آن که بیش‌ترین میزان 
تمرکز چاه‌ها و قنات‌ها را دارده از رسوبات آبرفتی تشکیل شده است که در مرکز و قسمت‌شرقی دشت بیش‌ترین 


میزان فرونشست دیده می‌شود و در سایر بخش‌ها به‌واسطه شیب بیش‌تر و نبود رسوبات آبرفتی زمینه برای وقوع 


TC‏ 9 بسیار کم‌تر از ساير بخش‌ها است. رسوبات جوان رودخانه‌ای و دشت با ۲۰ کیلومتر. مارن. 
توفهاء سنگآهک و مارن و ماسه‌سنگ توفی ۳۲ کیلومتر از مساحت پهنه بسیار يرخطر را شامل می‌شوند. 

بررسی معیار کاربری اراضی SE‏ از اين موضوع است که مناطق کشاورزی WAL‏ کیلومتر نواحی انسان‌ساز 
۹ کیلومتر دیم‌زار MIT‏ کیلومتر مناطق آبی ۹۸/۸ كيلومترء مراتع خوب ۱۵۵/0 کیلومتر» مراتع متوسط ۱۷۸/۸ 
شهرک صنعتی ۱۳۲/۹ کیلومتر مساحت مناطق بسیار پرخطر تا کم daz‏ را به خود اختصاص داده است. مناطق دارای 
کاربری اراضی کشاورزی آبی به دلیل بهره‌برداری از مناطق آب زیرزمینی» پتانسیل بالایی جهت وقوع فرونشست 
دارد. همچنین مناطق سکونت‌گاهی نيز به دليل ساخت‌وسازها و وجود ابنیه‌های بزرگ و سنگین مستعد وقوع 
فرونشست هستند. با توجه به ايتكه عمده محصولات کشاورزی زیر کشت دشت اردبیل گندم جو و سیب‌زمینی 
می‌باشد که نياز آبی بالایی برای رشد دارند و کاهش چشم‌گیر بارندگی و کمبود منابع آب‌های سطحی در سال‌های 
اخیر لزوم استفاده از منابع آب‌های زیرزمین را بالا برده و همین مساله سبب افت بیش‌تر منابع آب زیرزمینی در منطقه 


شده است. نتایج کاربری اراضی محدوده موردمطالعه و پراکنش فرونشست در اين پژوهش OUS‏ می‌دهد که به دلیل 
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رشد سریع جمعیت و توسعه صنعت و کشاورزی؛ تقاضای آب نیز بسیار زياد می‌باشد. بااین‌حال بهره‌برداری 
بیش‌ازحد آب‌های زیرزمینی به دلیل کمبود منابع آب‌های سطحی آشکار است که منجر به فرونشست زمین می‌شود 
و باعث آسیب به جادههاء يلهاء لوله‌های آب زیرزمینی و ساير تأسیسات شهری می‌شود. 


شکل ۷- تغییرات میانگین سطح ایستابی چاه‌های پیزومتری دشت اردبیل 


در دشت اردبیل قبل از دهه ٠١‏ و ۷۰ تعداد چاه‌های عمیق و نیمه‌عمیق در مقایسه با چاه‌های دستی يا نیمه عمیق 
انگشت‌شمار بود. درنتیجه برداشت از منابع آب زیرزمینی در حد محدود انجام می‌گرفته است. از سال ۷۸ به بعد 
استفاده از وسایل و دستگاه‌های حفاری, تعداد چاه‌های حفر شده در دشت اردبیل رو به افزايش نهاد و با رونق 
کشاورزی از دهه ۷۰ به بعد و به‌تبع آن افزایش بی‌رویه برداشت از سفره مذکور شروع به افت سطح آب زیرزمینی 
کرد. تداوم اين وضعیت در سال‌های بعدی به‌حصوص از سال ۱۳۸۵ به بعد باعث بحرانی شدن وضعیت اين دشت 
از بابت افت سطح آب و فرونشست شده است. کاهش آبدهی caola‏ چشمه‌ها و قنات‌ها و تغيير محل و کف‌شکنی 
چاه‌ها و درنتیجه تحمیل هزینه‌ای گزاف در اين خصوص به صاحبان چاه‌هاء افت کیفی منابع آب و شور شدن 
تدریجیی ایجاد محدودیت در مصارف کشاورزی» صنعت و شرب نشست‌زمین و کاهش ظرفیت حجم مخزن 
طبیعی از تبعات افت سطح آب‌های زیرزمینی در دشت اردبیل است. با استفاده از اطلاعات به دست آمده می‌توان 
كفت که قسمت شرقی و مرکزی دشت اردبیل بيش ترين تغیبرات را به لحاظ افت سطح آب زیرزمینی داشته‌اند كه به 
علت مساعد بودن شرايط کشاورزی و بالا بودن سطح زیر کشت محصولات با نیاز آبی بالا و تراکم چاه‌های عمیق 
حفر شده برداشت از آب زیرزمینی در منطقه ذکر شده فراوان بوده و به همین دلیل میزان افت سطح آب‌های 


زیرزمینی بالا می‌باشد. جهت بررسی وضعیت و تأثیر افت سطح آب‌های زیرزمینی و تأثیر آن بر میزان فرونشست 
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دشت اردبیل از اطلاعات ۳۹ چاه پیزومتری (اخذ شده از سازمان آب منطقه‌ای اردبیل)» طی سال‌های 51 تا ۹٩‏ (شکل 
۷ كه نتایج تحلیل از کاهش tA‏ درصدی سطح ایستابی را در سال ٩‏ نسبت به سال QUA TM‏ می‌دهد. بيشينه سطح 
ایستابی چاه‌های پیزومتری از ٤۵‏ متر به بیش از ۷۰ متر در سال ۹٩‏ رسیده است که حاکی از وخامت اوضاع آبخوان 
دشت اردبیل است. 

بررسی‌های زمینی مؤيد وجود تعداد کثیری شواهد و رخنمون‌های مورفولوژیکی ناشی از افت سطح آب‌های 
زیرزمینی و نشست زمین در دشت اردبیل می‌باشد. ازجمله می‌توان به درز و ترک‌های روی دیوار منازل (شکل ۸) و 
ترک و نشست جاده‌ها و خیابان‌های شهر اردبیل اشاره کرد و ترک‌هایی که در روی زمین‌های کشاورزی شکل گرفته 
سبب توسعه فرسایش شده است. به‌گونه‌ای که اين درز و ترک‌ها باعث کاهش کیفیت زمین‌های کشاورزی و از بين 


رفتن آن‌ها و آسیب به سیستم آبیاری شده است. 


i 


شکل ۸- ترک‌های شکل گرفته ناشی از فرونشست در دشت اردبیل 


غ- جمع‌بندی 

در اين پژوهش خطر فرونشست در دشت اردبیل» با استفاده از تلفیق روش‌های ANP‏ و MABAC‏ موردبررسی 
قرار گرفته است. نتایج مطالعه نشان colo‏ عوامل افت سطح آب (جزو عوامل هيدرو اقلیمی)» فاصله از آبراهه (جزو 
خوشه عوامل هيدرو اقلیمی» لیتولوژی (جزو خوشه زمین‌ساختی» به ترتيب با مقدار ۰/۲۲۱ ۰/۱۹۱ ۰/۱۵۲ 
پیشسترین ضريب وزنى را به خود احتصاص دادند. فرزین LS‏ و همکاران (۱۶۰۰) نيز در مطالعه‌ای که در آن به 
پهنه‌بندی خطر فرونشست در دشت جوین با استفاده از مدل تحلیل شبکه‌ای — فازی پرداخته‌اند؛ سطح جاه‌های 
پیزومتری» تراکم چاه‌ها و فاصله از رودخانه را به‌عنوان مهم‌ترین عوامل بیان کرده‌اند. شایان و همکاران (۱۳۹۵)» در 
مطالعه‌ای با عنوان پهنه‌بندی خطر فرونشست در حوضه‌ای آبخیز قره‌چای همدان به اين نتیجه رسیدند که تغییرات 
ناهمسان در داده‌های مورفومتری شبکه زهکشی مهم‌ترین فاکتور موثر در تشدید افت آب‌های زیرزمینی و وقوع 
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فرونشست در منطقه بوده و عوامل لیتولوژیکی و انسانی بهتنهايى نقش ندارند. با توجه به بررسی نقش و اهمیت هر 
يك از عوامل دخیل در شکل گیری فرونشست و نیز نتایج حاصل از بهكاركيرى روش تحقيق و تجزیه‌وتحلیل نهایی 
با استفاده از روش ماباک به‌عنوان یکی از روش‌های تحلیل تصمیم‌گیری چندمعیاره به ترتیب ۲6۵/۲۹ و ۳۷۰/۵۹ 
کیلومترمربع از محدوده مطالعاتی در طبقات بسیار پرخطر و پرخطر قرار دارد. با توجه به نقشة پهنه‌بندی فرونشست. 
دشت اردبیل دارای احتمال Jos‏ باه مناطق بسیار پرخط به‌طور عمده در شیب ۰ تا ۱۵ درصد. سازندهای نامقاوم و 
نیمه‌مقاوم. كه Gase‏ سازندهای سست کواترنری (رسوبات جوان رودخانه‌ای, رسوبات Sal‏ تراورتن. توف اهار 
مارن)» کاربری‌های زراعی: انسان‌ساخت و مراتع» محدوده چاه‌های پیزومتری و Grok‏ محدوده نزدیک گسل و در 
مناطق با فواصل دور از آبراهه‌ها قرار دارند. همچنین می‌توان كفت که نتایج حاصل از اين مطالعه. حاکی از توان GUL‏ 
UU ath c eme sse on rcs‏ وهای برعظر شین ad‏ در این ow‏ اراضی مهف که 
نباید کارهای غير اصولی از قبیل برداشت بی‌رویه آب‌های زیرزمینی. حفر چاه‌های عمیق و نیمه‌عمیق» کاشت 
محصولات با نیاز آبی بالا انجام گیرد. همچنین اين امر نیازمند شناسایی عوامل تأثي ركذار, اقدامات کاهش‌دهنده و 
اقدامات تطبیق دهنده ازجمله محدود يا ممنوعیت استخراج آب‌های زیرزمینی در مناطق در معرض فرونشست و 
اجرای مقررات سخت گیرانه جهت استخراج‌کنندگان منابع آب زیرزمینی در مناطق در معرض فرونشست و پایش 
پیوسته شاخص‌های فرونشست است. به‌طورکلی ایجاد چنین شرایط و وضعیتی در رویداد پدیده فرونشست و 
پیامدهای مخاطره‌آمیز و فاجعه‌بار آن برای كشورء قبل از هر چیزی عزم ملی برای پذیرش و مقابله با اين پدیده 
به‌عنوان یک معضل و مخاطره منتهی به يك فاجعه را طلب می‌نماید. موضوعی که درک أن از یک سو مستلزم 
افزايش مطالعات برای شناسایی کامل نواحی مستعد به‌صورت یک طرح ملی است. از طرفی لازم است زیرساخت 
لازم برای آگاهی‌رسانی و همکانی بودن اطلاع از نتایج و معضلات ناشی از اين مسئله برای مردم در جهت مشارکت 
آن‌ها در جهت جلب مشارکت مردم و مسئولان و تعدیل زیرساخت‌های ایجابی فراهم گردد؛ و نهايتاً برای پیشگیری 
و کاهش و تعدیل حجم دامنه فاجعه محتمل از اين مخاطره. دستگاه‌های ذی‌ربط تمامی مقررات و استانداردهای لازم 
ساخت‌وساز را برای نواحی در معرض فرونشست موردبازنگری قرارداده و اساس را بر شرایط نوين تعریف نماید. 
درنهایت می‌توان اذعان داشت با توجه به اين امر كه فرونشست زمین می‌تواند خسارت‌های جبران‌ناپذیر مالی و جانی 
به دنبال داشته باشد و به بسیاری از سازه‌های سطحی و زیرسطحی در مناطق شهری و حومة OI‏ در محدوده 
موردمطالعه سیب برساند. ازاین‌رو بایستی راهکارهای بلندمدت به‌صورت اصلاح روش‌های مدیریت منابع آب و 


راهکارهای کوتاه‌مدت به‌صورت. جلوگیری از برداشت سفره‌های آب زیرزمینی به اشکال مختلف و پایش مستمر 
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